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LONG-TERM  GOALS 

The  approach  that  describes  air-sea  exchange  as  interaction  between  a  flow  and  a  rough  surface,  now 
considered  traditional,  dates  back  to  the  1950s  (Charnock,  1956).  Within  that  approach,  the  complex 
exchange  processes  are  condensed  down  to  exchange  coefficients,  thus  delivering  computational 
efficiency  in  large-scale  numerical  modeling  of  air-sea  meteorology.  However,  inability  to  distinguish 
between  momentum  and  kinetic  energy  transferred  to  waves  from  those  transferred  to  currents,  as  well 
as  considerable  variance  in  the  experimental  estimates  of  the  drag  coefficients,  show  some  of  the 
applicability  limits  for  this  traditional  approach.  Clearly,  a  short-term  phase-resolved  wave  forecasting, 
a  goal  of  the  High-Resolution  Wave- Air-Sea  Interaction  project,  requires  a  more  detailed  mechanistic 
description  of  the  marine  boundary  layer  dynamics  with  a  special  focus  on  the  elements  distinctly 
introduced  by  the  compliant  interface  and  the  sea  surface  waves.  While  wave  dynamics  on  the  water 
side  has  already  been  reduced  to  a  computationally-intensive  numerical  problem  (Friehe  et  al.,  2007, 
section  III.B),  the  complexity  of  which  is  determined  by  the  number  of  nonlinearly  interacting  wave 
modes,  the  wind  driving  of  the  waves  on  the  other  hand,  is  less  understood.  Current  challenges 
include  gaps  in  theoretical  knowledge  and  in  techniques  for  numerical  modeling.  In  particular,  the 
observational  validation  for  most  of  the  wind-wave  interaction  mechanisms  proposed  so  far  in 
theoretical  works  is  lacking.  The  purpose  of  this  work  is  to  advance  our  understanding  on  these  open 
issues. 

OBJECTIVES 

A  goal  pursued  by  this  project  is  to  distill  both  field  observations  and  physical  analysis  into  a 
description  for  the  structure  and  dynamics  of  the  marine  atmospheric  boundary  layer  that  will  be 
suitable  to  incorporate  in  models  for  short-term  wave  prediction.  Considering  that  the  wave  breaking 
limits  the  waves  slope  and  thus  prevents  any  strong  nonlinearities  in  wind-wave  interaction,  as  a 
physical  framework  of  this  study  we  select  is  the  weakly  nonlinear  interaction  dynamics  proposed  by 
Hasselmann  (1965).  The  higher-order  mechanisms  in  that  hierarchy  have  been  a  subject  of  extensive 
theoretical  studies  in  the  last  four  decades,  yet  their  presence  and  intensity  has  never  been  addressed 
experimentally.  Some  non-expandable  mechanisms,  such  as  flow  separation,  have  received  renewed 
interest  in  light  of  the  work  of  P.  Sullivan,  NCAR.  Detection  of  such  mechanisms  in  measurements 
will  be  among  our  important  objectives. 

From  applied  perspective,  we  focus  on  experimental  determination  of  the  wind  input  to  the  waves,  an 
essential  physical  factor  in  waves  evolution.  Our  analysis  has  lead  us  to  conclude  that  a  widely  adopted 
approach  of  exponential  extrapolation  to  obtain  atmospheric  surface  pressure  from  pressure 
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measurements  in  the  wind,  has  a  clear  tendency  of  overestimating  the  surface  pressure  and  may  have 
caused  substantial  uncertainty  in  previous  estimates  of  wind  input  to  the  waves.  Here  we  propose  an 
alternative  that  is  free  of  the  deficiencies  of  exponential  extrapolation. 

APPROACH 

Meteorological  and  wave  data  acquired  in  June  2010  over  the  open  ocean  will  be  analyzed,  along  with 
data  from  previous  experiments,  to  assemble  a  complete  description  of  the  marine  boundary  layer 
dynamics.  As  the  2010  data  cover  mostly  conditions  of  high  winds  and  developed  wave  fields,  data 
from  previous  experiments  will  complement  our  study  with  information  on  low  and  moderate  winds. 
The  wind  input  results  will  be  used  to  calibrate  numerical  models  of  wave  evolution  and  marine 
boundary  layer  dynamics.  The  atmospheric  pressure  measurements  will  be  extrapolated  to  the  surface 
in  a  novel  and  more  rigorous  way,  to  evaluate  the  direct  energy  input  to  the  waves  occurring  through 
form  drag.  The  wind  velocity  measurements  will  be  analyzed  to  identify  wave-turbulence  interaction 
and  its  contribution  to  the  wave  growth.  Recent  numerical  results  reported  by  Peter  Sullivan,  NCAR, 
have  predicted  instances  of  air  flow  separation  in  cases  of  extreme  wave  slope  and  surface  roughness. 
We  will  seek  to  recognize  signs  of  such  phenomena  in  the  data  and  estimate  their  influence  on  the 
wave  field  evolution. 

The  PI  closely  collaborates  with  other  members  of  the  project's  team.  The  experimental  component  is 
carried  out  with  Carl  Friehe  and  Jesus  Ruiz-Plancarte  (UCI).  The  work  on  data  assimilation,  surface 
and  boundary  layer  air  flow  modeling  will  be  done  with  Eric  Terril  (UCSD)  and  Peter  Sullivan 
(NCAR). 

WORK  COMPLETED 

An  extensive  preparation  for  the  field  experiment  included  designing  and  building  deployment  gear, 
testing  of  instruments  and  developing  of  data  acquisition  and  data  processing  algorithms  and  software. 
Starting  June  1,  2010  FLIP  moved  to  the  site  selected  for  the  experiment  Northwest  of  San  Francisco, 
at  33.3365N  and  123.4289W.  Site’s  climatology  of  high  winds  and  large  waves  significantly  extended 
the  time  necessary  for  deployment.  Atmospheric  conditions  and  sea  state  included  wind  gusts  up  to  2 1 
m/s  and  significant  wave  heights  up  to  18  feet.  The  lower  part  of  the  mast  has  been  repeatedly 
submerged  with  some  repairable  damage.  Among  the  essential  components  of  the  meteorological 
array  were  the  instruments  measuring  airflow  velocity  (sonic  anemometers,  cups,  vanes),  pressure 
sensors,  and  the  units  for  GPS  tracking  and  inertial  navigation.  Total  of  57  units  were  deployed  and 
142  signals  were  acquired  continuously,  most  of  them  listed  in  Table  1.  The  instrumented  mast  and  the 
persistent  whitecaping  on  the  sea  surface  are  shown  in  Figure  1.  We  gratefully  acknowledge  Ken 
Melville,  UCSD,  and  Fabrice  Veron,  Univ.  of  Delaware,  who  kindly  shared  the  signals  from  two  laser 
altimeters  registering  wave  elevation.  Strong  current  first  led  to  excessive  drag  and  tilt  on  FLIP  that 
later  caused  the  main  mooring  line  to  snap.  Measurements  had  to  be  discontinued  on  June  19,  before 
the  scheduled  end  date,  followed  by  accelerated  redeployment. 

RESULTS 

A.  Wave  influences  in  the  pressure  and  velocity  fields. 

Out  preliminary  analysis  of  the  data  set  was  in  part  inspired  by  the  recent  findings  of  Peter  Sullivan, 
NCAR,  whose  modeling  results  indicated  a  significant  distinction  between  the  structures  of  velocity 
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and  pressure  fields  in  the  wavy  boundary  layer  (Sullivan  et  al.,  2010).  Namely,  modeling  results  find 
that  the  wave  influence  in  pressure  is  more  pronounced  than  in  velocity,  a  fact  with  potentially 
significant  implications  on  the  dynamics  of  wind-wave  coupling.  Such  findings  agree  with  earlier  field 
observations,  e.g.  from  the  RED  and  CBLAST  experiments.  A  consistent  interpretation  is  pursued 
within  this  project. 


Instruments 

Quantity  measured 

Height/Location 

12  Sonic  Anemometers 

Wind  velocity, 

1 0  on  the  mast,  2  on  the  boom 

10  RMY  Prop.  Anemometers 

Air  temperature 

Wind  speed  &  direction 

At  1 0  levels  on  the  mast 

1 0  Pressure  Instruments 

Atmospheric  pressure 

8  on  the  mast,  2  on  the  boom 

2  Wave  Wires 

Sea  surface  elevation 

At  the  mast  and  at  the  middle  of  the 
boom 

2  Laser  Altimeters 

Sea  surface  elevation 

At  the  mast  and  at  the  middle  of  the 
boom 

Krypton  Hygrometer 

Atm.  humidity 

On  the  mast,  at  1 1.3m  height 

Infrared  SST 

Sea  surface  temperature 

Leosphere  Windcube 

Wind  Profiles  in  the  ABL 

FLIP’S  Gyroscope 

Boeing  CMIGITS  III 

Heading 

Platform  Motion 

On  the  mast,  at  1 8m  height 

GPS  Unit 

GPS  location  and  timing 

8  Thermistors 

Air  temperature 

5  on  the  mast,  3  in  the  water 

Oxford  Inertial  Nav.  Unit 

GPS  position,  Motion 

On  the  boom 

Table  1.  List  of  instruments  deployed  during  the  field  experiment  in  June  2010.  The  Laser  Altimeters 
were  kindly  shared  by  Ken  Melville  (UCSD)  and  Fabrice  Veron  (Unix.  Delaware). 
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Figure  1.  The  instrumented  mast  on  the  port  boom  (left) 
and  the  persistent  whitecaping  on  the  sea  surface  at  high  winds  (right). 


One  could  observe  that  for  uncorrelated  waves  and  turbulence  the  correlation  function  of  the  measured 
velocity  B"'  (q ,  r2 )  and  pressure  Bp  (f ,  r2 )  split  into  turbulent  Bfr)  ,B'p(r)  (where  r  =  rx  -r2)  and 

wave-induced  terms  B”(rx,r2)  and  Bwp(rx,r2),  i.e.  Bnfp]  (7]  ,r2)  =  B[vp]  (r)  +  Bfp,  (q ,r2) .  Recalling  that 
the  structure  D{r )  and  correlation  B{r )  functions  of  the  same  random  field  are  related  as 


B(r )  =  f?(0)  -  ^  D(r  ) ,  and  that  for  the  inertial  subrange  of  developed  turbulence  the  turbulent  structure 

functions  scale  as  Bfr)  =  Cv2r2/3  for  velocity  and  B'p(r)  =  C2pr4,i  for  pressure,  one  arrives  to  these 
forms  of  the  correlation  functions 

B:(fltr2)  = 

B:(rl,r2)  = 


B‘M-'-Cyn  +B;yr i) 


1 


B'A  0)--c> 


2 _ 4/3 


+  B:(fl3r2) 


These  suggest  that  the  turbulence  in  pressure  decays  faster  with  separation  distance  (exponent  (4/  3) ) 
Than  it  does  in  velocity  (exponent  (2/3) ).  Consequently,  the  influence  of  the  boundary,  i.e.  the  waves, 
will  be  more  pronounced  and  propagate  further  in  the  flow  in  pressure  than  in  velocity,  as  has  been 
observed  in  low  to  moderate  wind  conditions.  At  high  winds,  however,  two  circumstances  occur 
simultaneously:  (i)  there  is  a  lowering  of  the  critical  layer  in  the  air  flow,  which  roughly  determines  the 
vertical  extent  of  the  wave-induced  pressure  and  (ii)  there  is  increased  turbulent  intensity. 
Consequently,  turbulence  also  starts  to  dominate  the  pressure  signal  and  no  clear  wave  signature  is 
observed,  as  indicated  by  the  pressure  spectrum  in  Figure  2. 
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B.  Determining  wind-to-wave  energy  input  from  measurements  of pressure  in  the  airflow 

Pressure  distribution  on  the  sea  surface  detennines  the  rate  of  wind-to-wave  energy  transfer  (f  f'j . 

However,  direct  measurements  of  the  atmospheric  pressure  on  the  surface  are  impractical,  considering 
that  the  sensitive  instrument  will  have  to  be  protected  from  wetting,  and  that  the  film  serving  that 
purpose  will  distort  the  pressure  readings.  Also,  avoiding  spurious  pressure  readings  from  dynamic 
pressure  requires  a  specific  pressure  port  (inlet)  whose  disks  should  be  kept  parallel  to  the  mean  flow 
(Hristov,  2008).  This  requirement  is  difficult  to  meet  when  both  the  position  and  the  orientation  of  the 
instrument  are  vigorously  driven  by  the  surface's  motion.  Furthermore,  as  the  pressure  on  the  surface 
does  not  directly  detennine  the  structure  of  the  airflow  above,  information  on  the  surface  pressure 
carries  little  or  no  information  on  the  wind-wave  generation  dynamics.  Instead,  pressure  measurements 
in  the  airflow  are  free  of  the  experimental  limitations  mentioned  above  and  can  reveal  more  of  the 
wind-wave  coupling  mechanism,  yet  for  estimating  wind-wave  energy  transfer  rate  these  have  to  be 
extrapolated  down  to  the  air-water  interface.  The  extrapolation  method,  essential  for  proper 
determination  of  the  wind  input  to  the  waves,  is  proposed  for  this  project  and  is  outlined  below. 

The  velocity  and  pressure  variation  induced  in  the  water  by  deep-water  waves  exhibits  vertical 
exponential  decay.  In  previous  works  (Donelan,  1999),  (Hasselmann  and  Bosenberg,  1991)  authors 
have  assumed  an  analogous  behavior  for  the  wave-induced  motion  on  the  air  side  and  have  used  that 
assumption  to  extrapolate  pressure  measurements  within  the  airflow  down  to  the  interface.  This 
assumption,  however,  is  inconsistent  with  the  presence  of  singular  point  in  the  air  flow  and  may 
considerably  overestimate  the  wind-wave  energy  input.  Namely,  the  height  where  the  mean  wind 
speed  matches  the  phase  speed  of  the  wave  carries  dynamic  significance  and  the  vertical  structure  of 
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Figure  2.  A  spectrum  of  the  atmospheric  pressure  (upper  plot)  and  a  wave  spectra  produced  by  a 
wave  wire  (blue)  and  laser  altimeter  (green),  shared  by  Ken  Melville  (UCSD)  and  Fabrice  Veron 
( Univ .  Delaware)  (lower  plot).  At  the  high  wind  conditions  encountered  throughout  the  experiment, 
the  wave  signature  in  the  pressure  is  weak,  a  consequence  of  the  lowered  critical  layer  and 

increased  turbulent  intensity. 


the  wavy  boundary  layer  is  profoundly  influenced  by  that  critical  layer  singularity.  The  critical  layer 
theory  (Miles,  1957),  (Hristov  et  al.,  2003)  properly  accounts  for  the  role  of  the  matching  height  and 
shows  vertical  behavior  for  the  pressure  that  is  distinctly  different  from  exponential.  Specifically,  for  a 
range  of  wave  ages,  the  pressure  shows  virtually  no  vertical  decay  and  change  in  phase  up  to  the 
matched  height.  Within  the  critical  layer  theory,  the  stream  function  </){z)  of  the  wave-induced  flow 


satisfies  the  Rayleigh  equation  £,{(/)"  -(f)-  f'(f>  =  0 ,  with C  =  log(z / zc) .  Expressing  the  pressure 

dtj) 


r  j±\ 


through  wave-age  parameter  dfl  =  -ck  /  u,  and  the  solution  <f>(z) ,  i.e.  p  =  (p„g)  - 


Q 


4 

v  v 


dg 


-<t> 


V 


where  Q  is  the  Charnock’s  roughness  parameter,  one  could  relate  the  wave-coherent  pressure  at  a 
given  height  p(g)  to  the  pressure  at  the  surface  p(f  t)  as  p(fi})  =  T(f,  f  ^pif) .  The  transfer  function 
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r(££0)  takes  the  form: 


/ 

\ 

-  1 

£0  , 

<?  =  -fo  y 

fo-f 


f  d<p^ 
v  v 


d % 


-  4> 


Application  of  this  transfer  function  to  sampled  time  series  requires  the  discrete -time  filter  counterpart 

+  b,X 


b  b  z~ ^ 

of  that  function,  commonly  in  z  -transfonn  representation  F(z)  =  0  1 


1  +  atz  H — amz~ 


.  Such  a 


counterpart  filter  can  be  obtained  by  selecting  the  order  of  the  filter  (n,  m )  and  requiring  that  the 
discrepancy  s  between  the  amplitude  and  phase  responses  of  the  complex-valued  analog  prototype 
7'(f;f0)  and  the  discrete  form  F(z) 


s{bQ,bx,a0)=  Yj  \(l  +  a^y (a)  -  (b{  +  b^J 
n 

is  minimized.  Choosing,  for  instance,  n  =  1,  m  =  1 ,  and  necessitating  that  the  derivatives  at  the 
minimum  should  vanish,  leads  to  these  equations  for  detennining  the  coefficients  \bQ,b{,a{}: 

ff  =  Z  K  +  f "  +  <’""X  -  (e"r  +  -  {t  +  r*)]=  o 

db0  n 


|f  =  £  [(«'"  +  «■'"  >0  -  (t  +  T‘  k  -  (e"T  •  +  <r'*T)]  =  0 


—  =  X  [(e'wr  +  e~i(OT*)b0  -  (t  +  T*\  +  2  TT  *  a,  +  (e*®  +  e-'®  )zT  *]  =  0 

Above,  <y  is  the  frequency  measured  in  radians  per  sample.  The  set  of  coefficients  {b(>  ,b],a] }  depends 
on  the  distance  of  the  instrument  from  the  surface,  the  wind’s  friction  velocit  u„ ,  and  the  pressure 
signal’s  sampling  frequency.  Considering  that  the  pressure  pm  measured  in  the  airflow  broadly  consists 
of  a  turbulent  pt  and  a  wave-induced  p  terms,  i.e.  pm  =  p,+p ,  the  filter  F{ico)  applied  to  the  wave- 
induced  component  p(^)  produces  the  wave-coherent  pressure  at  the  surface  p(gtl ) ,  suitable  for 
calculating  the  wind-to-wave  energy  transfer  rate  (p(<^0)rj)  ■  For  computational  efficiency,  one  may 
skip  the  separation  pm  into  a  turbulent  pt  and  a  wave-induced  p  parts.  Using  either F(ico)  0  pm(Jg)  or 
p(40)  should  produce  the  same  energy  transfer  estimate,  i.e.  (pig,,  )>])  =  ( rjF(ia> )  0  pm (g)) .  Since  the 
turbulent  component  p,(4)  and  its  (unphysical)  filtered  counterpart  F{ico)  0  p,(g)  are  both 
uncorrelated  with  the  waves,  i.e.  ( F(ia> )  0  pt  (£),  77)  =  0  ,  F(ico)  0  p,  {Jg)  will  not  contribute  to  the 
energy  transfer  rate. 


IMPACT/APPLICATIONS 


The  results  of  this  research  are  expected  to  advance  the  basic  science  of  the  air-sea  interaction  and  will 
be  applied  to  operational  models  for  short-term  wave  modeling  and  forecasting.  The  information  on 
the  structure  and  dynamics  of  the  marine  atmospheric  boundary  layer  (MABL)  and  the  statistics  of  the 
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ocean  surface  will  advance  the  description  and  modeling  of  signal  propagation  over  the  ocean.  The 
profound  physical  similarities  between  propagation  of  radar  signals  over  the  ocean  and  acoustic  signals 
in  the  water  will  extend  possible  applications  to  the  acoustic  domain. 


RELATED  PROJECTS 

The  PI  is  unaware  of  any  related  projects. 
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